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Abstract 

We investigate the general case of the photon distribution of a two-mode squeezed vacuum 
and show that the distribution of photons among the two modes depends on four param- 
eters: two squeezing parameters, the relative phase between the two oscillators and their 
spatial orientation. The distribution of the total number of photons depends only on the 
two squeezing parameters. We derive analytical expressions and present pictures for both 
distributions. 


1 Introduction 

Squeezing the quantum fluctuations of the radiation field has been demonstrated experimentally 
using various optical systems [1]. Most of them rely on two-mode squeezing [2, 3, 4]. Therefore 
the properties of a two-mode squeezed state have been studied extensively [5, 6]. However it 
appears, that the photon statistics of such a state has not been investigated in all details. Only 
some particular cases have been considered (see for example Refs. [/] and [8]). In the present 
paper we therefore extend these considerations to an arbitrary two-mode squeezed vacuum and 
address the following questions: (i) What is the most general case of the squeezed vacuum of a 
two-mode oscillator, and how many independent parameters are needed to describe this state? 
(ii) What is the photon statistics in this state? This is of interest in the context of the degenerate 
parametric amplifier [9] since the mathematical structure of the two-mode distribution function 
coincides with the transition probability function of a two-dimensional parametric oscillator [10]. 
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2 general case of the squeezed vacuum of a two- mode 
oscillator 

We start our considerations with the question (i). The linear canonical transformation 

(a',a;b\b)^{A\A;B\B) (1) 

of the creation operators a*, 6* and the annihilation operators a, b to their new counterparts A\ 
B t and A, B corresponding to the two modes suggests a set of 10 parameters. Indeed, the generic 
transformation (1 ) of the four-dimensional vector (a*, a, b\ b) into (A\ A. B\ B ) via a 4 x 4 matrix 
brings in 16 complex values, that is 32 parameters. Since a and a\ b and b * ,A and A\ B and B * 
are hermitian conjugate of each other, only half of these parameters are independent. Therefore 
the number of parameters reduces to 16. The condition to preserve the commutation relation 
[A, B] = 0, [A, /?*] = 0, [y4, A^] = 1, and [B,B^] = 1 provides additional constraints. Since the 
commutator of an operator with its hermitian conjugate is always real, the last two conditions 
provide only two constraints, whereas the first two conditions must hold for the real and the 
imaginary part separately. This decreases the number of parameters by six leaving us indeed with 
10 parameters. 

This finding is in accordance with the results of group theory. Reference [11] shows, that 
the group of rotations and squeezing of the four-dimensional phase space, that preserves the 
phase space volume of the two degrees of freedom, consists of 10 generators. Hence there are 10 
parameters determining the elements of the real symplectic group Sp(4,R). 

But what is the physical meaning of these parameters? They are associated with rotation and 
squeezing transformation of phase space of two oscillators. Generic rotations of four-dimensional 
phase space are described by 6 parameters but the symplectic rotations in four-dimensional space, 
that is the transformation preserving the commutation relations, are described by 4 parameters. 
They can be represented as a sequence of four rotations. The first rotation given by the trans- 
formation where M — i(a^b — b^a) is the angular momentum operator, corresponds to the 

rotation in the coordinate space by the angle fa. The second and the third rotations are given by 
the operator e*^ ata and e ,xfct6 , and correspond to rotations by the angles ip and \ in the respec- 
tive phase-spaces of the two oscillators. The last rotation again can be taken in the form e tA ^ 2 . 
Thus the angles <f >\ , <f) 2 , ip and \ are chosen to be the parameters of the symplectic rotation. 
The general symplectic transformation also includes squeezing. We can represent the generic real 
symplectic transformation (1) as consisting of three consecutive transformations: symplectic rota- 
tion (4 parameters), followed by independent squeezing of the two modes given by the squeezing 
operators 

5j = e n(° 2 -at 2 )/2 5 5 2 _ e Mb 2 -b< 2 )/2, (2) 

(2 squeezing parameters ?q and r 2 ) and followed by another real symplectic rotation (4 parameters 
more). But how many of these parameters govern the two-dimensional squeezed vacuum state? 
The first four parametric rotation acting on the completely symmetric vacuum state leaves this 
state unchanged. Hence 2 squeezing parameters and 4 parameters associated with the second 
rotations define the two-dimensional squeezed vacuum state. In addition we can include an overall 
quantum phase factor t lp of this state. An explicit calculation of the general case of a two-mode 
squeezed vacuum wave function is given in Ref. [12]. 
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3 Photon statistics for the total number of photons 


We now address the question: how many of these parameters govern the photon statistics ol 
such a two-dimensional squeezed vacuum state? In particular, how many of them determine the 
probability of (i) counting a total number n of photons in the two modes, and (ii) counting the 
number iii and n 2 of photons in the individual modes? 

In this section we consider the case of (i) and ask for the total number n of photons in both 
modes, which corresponds to the surface of a four-dimensional sphere, + x\ + p\ + orij) = n, 
centered at the origin of the four-dimensional phase space. The rotation of the phase space does 
not alter this sphere and hence only two squeezing parameters are essential. The probability of 
counting the total number n of photons for the case of two independently squeezed oscillators 
reads 

w n (s u s 2 ) = it Wn n,(s 2 ), (3) 

Til =0 

where W nj (sj) is the one-dimensional photon statistics [13] 


W^) - 


0 for nj odd 

VI - ~i S 7' 2 2 ~ UJ (n"/ 2 ) f ° r n i eVen 


of the squeezed vacuum wave function 


where 

Equation (4) reduces Eq. (3) to 


,r ,/2 


'M*j) = ^TTT ex P( — 2 e rj ' T i) 


Sj = tanli 2 (?•_,). 


(4) 

(5) 

( 6 ) 


(^1 ■) ^2 ) 


w 2t (s u s2) = 

^2fc+l(-Si, S 2 ) = 0 


(7) 


where the normalization factor Af reads 


Af — y/\ — S\ \J\ — S 2 - ( 8 ) 

The odd terms of W n vanish because of Eq. (4) . The sum in Eq. (7) has been calculated in [12] 
and the probability of counting n = 2k photons then reads 

W 2k (s u s 2 ) = Afs k 2 2 Fi{-k, 1/2, 1; 1 - s,/s 2 ) (9) 

where 2 F\ describes the hypergeometric function. For the special case of identical squeezing in 
the two modes, that is, ,s 1 = .$ 2 = .s, Eq. (9) yields 

W 2k {s) = .Vs k . ( 10 ) 
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In Fig. 1 we show the photon statistics (9) for various magnitudes of the squeezing parameters 
S\ and $2- The solid and dashed curves correspond to weak and strong symmetric squeezing 
($i = $ 2 ), respectively. In accordance with Eq. (10) the photon distribution then displays an 
exponential dependence. Stronger squeezing results in a higher amount of quanta involved. The 
dotted curve shows the photon statistic for an asymmetric squeezing. 



Fig. 1. Probability of counting n = 2k photons in a two-mode squeezed vacuum 
as given by Eq. (9). For symmetric squeezing (dashed line for s\ = s 2 = 0.5 and solid 
line for S\ = s 2 = 0.99) the curve is a straight line and hence exponential whereas for 
asymmetric squeezing (dotted line for S\ = 0.5 and s 2 — 0.99 ) the photon statistics 
is non-exponential. Here we have not specified the distribution of the n = 2k photons 
among the two modes. Note that W n = 2 k + 1 = 0 which we have omitted for simplicity. 


4 Photon statistics in the individual modes 


We now turn to the second case and calculate the photon statistics W (n lt n 2 ) = | < Uj, UjI'I',, > | 2 
in the individual modes. For this purpose we start from the generic expression of a squeezed 
vacuum state |'I' ag >, which we produce in the following way: 


i,x 2 ) 


= e' fi e~'* Kl f( 7 ) Si(r,) S 2 (r 2 ) * 0 , 0 (x u x 2 ) 





(ii) 


with the vacuum wave function of the two-dimensional oscillator ^o,o(* r iT r 2) = ^0(^1) ^o( r 2)i the 
angular momentum operator M = i (a^b — tfa) and the operator f(q') = exp — a^a)) which 
describes a mutual pha^e shift 27 between the oscillators. The wave function of a two-dimensional 
squeezed vacuum state is a Gaussian described by the three complex numbers /t, B . and C, which 
are functions of the four real parameters 77, r 2 , and 7. Their explicit dependence is given in 
[12]. One also finds there the generic two- mode squeezed vacuum wave function depending on two 
more parameters that do not affect the photon distribution. 
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We contract Eq. (11) with the probability amplitude of the photon energy states and 

'I’nii-Vi) and arrive after calculating the resulting double integral [12] at 


H 'f r r \ - 8 exp(-|ln(n,!/n 2 !)|) r 2x,/ 2 .. 

l!(?,, ’ Wa) ~ i(2^ + l)(2fi + l)-4C 2 | Ml5, “ Cl) X 


\(2A + 1)(25 + 1)-4C 2 | 

AAB + 2A-2B-1- 4 C 2 
4AB - 2A + 2B -1 - 4 C 2 

4A5 - 2/1 - 25 + 1 - 4C 2 

4,45 + 2.4 + 25+1 - 4C 2 

/ 

i|»l — Tl2|/2 


(»1 “«2)/2 


( n !+« 2 )/2 


X 


-4C 


( 12 ) 


|ni+U2l/2 W4A5 + 2+ + 25 + 1 - 4CV4C 2 - 4.45 + 2,4 + 25-1 


Here P k denotes the associated Legendre Polynomial. As in Eq. (7) the total number of photons 
n = 1 ) \ +?} 2 must be an even number hence 7i\ and n 2 are either both even or both odd. Otherwise 
the photon distribution function vanishes. 

In Fig. 2 we display the probability to find n 1 photons in mode 1 and n 2 photons in mode 2. 
This is the generic case of the photon distribution of a two-mode squeezed vacuum. It depends 
on four parameters: two squeezing parameters, the orientation of the distribution function with 
respect to our laboratory system, and the correlation between the two modes. Beside the even- 
odd oscillations the maxima lie on curves which are symmetric with respect to the main diagonal 
7?i = ii‘ 2 . This behavior is similar to the distribution function of a displaced two-mode squeezed 
state discussed in Ref. [7, 8]. 



Fig. 2. Probability of counting and n 2 photons in the two-mode squeezed 
vacuum calculated via Eq. (12). We have chosen 77 = 3, = 5. 6 = f . and 7 = — . 

The wavy structure of this distribution is confined to an angle in the photon number 
plane. 
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5 Conclusion 


We conclude by summarizing our main results. The wave function of the two-mode squeezed 
vacuum depends on 6 parameters (besides the phase factor). However only 4 of them, - the 
squeezing parameters r } and 7-2 of the two modes, the phase difference 7 between the two oscillators, 
and the rotation of the reference system by the angle <j) - manifest themselves in the generic case in 
the distribution of the photons among the two modes. This distribution function can be expressed 
explicitly in terms of Legendre polynomials. Only two parameters 77 and r 2 govern the distribution 
of the total number of photons, which we express explicitly in terms of a hypergeometric function. 
In conclusion we want to make the remark that similar consideration for iV-mode squeezed vacuum 
state shows that the photon distribution W(tii, n 2 , ...,nyv) depends on N 2 parameters. 
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